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Abstract 

We study the efficiency at maximum power, r]m, of irreversible quantum Carnot engines (QCEs) 
that perform finite-time cycles between a hot and a cold reservoir at temperatures and Tc, 
respectively. For QCEs in the reversible limit (long cycle period, zero dissipation), rjm becomes 
identical to Carnot efficiency rj^ = 1 — For QCE cycles in which nonadiabatic dissipation and 
time spent on two adiabats are included, the efficiency r/^ at maximum power output is bounded 
from above by and from below by In the case of symmetric dissipation, the Curzon- 

Ahlborn efficiency r/^^ = 1 — is recovered under the condition that the time allocation between 
the adiabats and the contact time with the reservoir satisfy a certain relation. 
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I. INTRODUCTION 



The concept of Carnot efficiency is of paramount importance in thermodynamics, since 
the Carnot cycle is the most efficient heat engine cycle allowed by physical laws. When 
the thermodynamic second law states that not all the supplied heat is applied to producing 
work, the Carnot efficiency presents the limiting value on the fraction of the heat which 
can be so used. Although the quasistatic Carnot cycle has the highest efficiency, it outputs 
zero power because it takes infinite time to output a finite amount of work. By contrast, 
Curzon and Ahlborn l| considered a finite-time Carnot cycle under the assumption of 
endoreversibility that irreversible processes occur only through these heat exchanges, they 
obtained the efficiency r)^^ at maximum power output as 



CA 



(1) 



where Th and are the temperatures of the hot and cold heat reservoirs, respectively. 
The Curzon- Ahlborn (CA) paper has triggered the development of research into finite time 
thermodynamics js iQ]. The above ri^ is usually called the CA efficiency, describing the 



efficiency of several engine models 
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1 120[ and of actual thermal plants 
very well. The CA efficiency has been found to be a universal result in the case of the low, 
asymmetric dissipation, by optimizing power output with respect to time allocation when 
time durations in adiabats and nonadiabatic phenomenon were ignored 
Great efforts have been devoted to the study of quantum heat engines 
beginning with the concept of quantum heat engine introduced by Scovil and Schulz-DuBois 
30| . Quantum heat engines differ from classical counterparts mainly in the following three 
respects: fi) the working substance is composed of quantum matter such as spin systems 
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261], harmonic oscillator systems 
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2lH24|. cavity quantum electrodynamics systems [2J, |27H29||. etc. (ii) The state of the 



12 L two- level or multilevel systems 



12| . (iii) Quantum equations 



system is depicted by a quantum-mechanical operator, and the thermodynamic observables 
are associated with the expectation of values of operators [lO 
of motion are used to describe the time evolution of the observab 
which can avoid using phenomenological heat transfer laws 
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es in quantum heat engines. 
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The previous literature discussed the heat engine models in the sudden limit in which 
the adiabatic process is a spontaneous switching and thus the time allocation on adiabats 
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is negligible [8|, |16|, |17[. Thus it is significant to study more general models in which the 
"adiabatic" process ( we take the two corresponding processes as two quantum "adiabats" 
throughout the paper) takes finite time as well become nonadiabatic 3l| . During a quantum 
adiabatic process, the variation of the eigenspectrum (quantum state) of the system must be 
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32| can apply. Otherwise, nonadiabatic 



so slow that the quantum adiabatic theorem 

dissipation (e.g., inner friction 26|) occurs because of rapid change in the energy level 

structure of the quantum system. Particularly nonadiabatic dissipation has been found to 
have a profound infiuence on the performance of quantum heat engines [ll[[l2[[26|. Including 
nonadiabatic dissipation is therefore essential for more realistic models of quantum heat 
engines. 

In this paper, we study the efficiency at maximum power output of QCEs performing 
finite time cycles, in which the time of any adiabat and nonadiabatic dissipation are consid- 
ered. We assume that the external parameter affecting the energy spectrum varies at a small 
but fixed speed which, however, may not be slow enough and thus to cause nonadiabatic 
phenomenon. We derive the cycle period that consists of times spent both on the two quan- 
tum isotherms and on the two quantum adiabats. We show that the efficiency at maximum 
power output converges to an upper and a lower bound in the limits of extremely asymmetric 
dissipation. Based on the low- dissipation assumption that the irreversible entropy produc- 
tion in a thermodynamic process is inversely proportional to the time required to complete 
that process, our approach similar to that of the classical thermodynamics predicts that the 
CA efficiency turns out to be an exact and universal property for QCEs operating under the 
conditions that the dissipation is symmetric and the time allocation between the adiabats 
and the contact time with the reservoir satisfy a certain relation. 



II. EFFICIENCY AT MAXIMUM POWER OUTPUT 

We consider a quantum system whose Schrodinger's equation is given by H\un) = En\un)-i 
where if, and En are the Hamiltonian of the system, its nth eigenstate and eigenenergy, 
respectively. The internal energy U reads U = EnPn, where P„ is the mean occupation 
probability of the nth eigenstate and obeys the canonical distribution P„, = Ig^^n/kBT 
equilibrium, with the canonical partition function Z = ^-En/ksT^ Derivation of U leads 
to the first quantum thermodynamic law dU = J^n^ndPn + J^n^ndEn- Analogous to the 
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classical thermodynamic first law, the first law of thermodynamics in quantum-mechanical 



systems is |23|-|25| dU = dQ +dW, in which (fQ = J^n^ndPn a.nd dW = J^n^ndEn depict 
the heat exchange and work done, respectively, during a thermodynamic process. Note 
that EndPn is associated with the heat exchange because dQ = TdS with the entropy 
S = —ksPn^T^ Pn- Motivated by the definition of the generalized force F for a classical 
system, we define analogously the force for a quantum system as F = J2n ^^qx^^ ? where 
X is the external parameter (generalized coordinate corresponding to the force F) [23 
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33| . Here the force F and generalized coordinate X are state variables 
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quantum versions of the classical pressure Pr and volume V, respectively. 
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FIG. 1: Schematic diagram of an irreversible QCE cycle in the plane of the external parameter 
X and force F{X) (a) and of the Temperature T and entropy S (b). The values of the external 
parameter X and of the entropy S at the four special instants are indicated. 



The generalized force F{X) alters the generalized coordinate X that affects the eigen- 
spectrum of the system during a thermodynamic process. The quantum Carnot cycle 
1— i'2— 7-3— i'4— T-lis drawn in the {F, X) plane, which is the quantum version of 
the classical {Pr, V) plane (See. Fig. [1^). During two quantum isothermal processes 1 — )■ 2 
and 3 — )■ 4, the working substance is coupled to a hot and a cold heat reservoir at constant 
temperatures Th and Tc, respectively. We apply dQ = TdS directly to the calculation of the 
heat exchange dQ in any quantum isothermal process. Let S{Xi) and Xi be the entropies 
and the external parameters at the instants i with i = 1,2,3,4, the heat amount Qh ab- 
sorbed from the hot reservoir and the heat amount Qc released to the cold reservoir are, 
respectively, Qh = Th[S{X2) - S{Xi)] and Q^ = T^\[S{X4) - ^(Xg)]]. When the Carnot 



cycle is reversible, the system couples to the heat reservoir for a sufficiently long time until 
the system remains infinitesimally close to equilibrium all along the cycle, and in order for 
the adiabatic theory to remain valid, the time scale of the ch ang e of the quantum state must 



be much larger than that of the dynamical one, ^ E/h 
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32|. For the reversible cycle 



where S{X2) = S{X^) and S{Xi) = Si^X^), we recover the Carnot efficiency r^^ = 1 — 
which is independent of the properties of the quantum working substance. We denote by 
ti2 (^34) the time durations during which the system is coupled to the hot (cold) reservoir 
along a cycle. In the branch 2 — )■ 3 (4 — 1), the working substance is decoupled from the 
hot (cold) reservoir, and the entropy changes from 5*2 to S3 (5*4 to 5*1) during a period ^23 ( 
^41), as shown in Fig. [T)d, where Si = S{Xi) with i = 1,2, 3, 4. 

Let us consider a QCE under finite-time operation. Finite-time cycles move the system 
away from the equilibrium, leading to irreversibility of the engine. Although the system 
needs no close to equilibrium during the isothermal process, the system remains in an equi- 
librium state with the heat reservoir at special instants i with i = 1,2,3,4. Under such a 
circumstance, the thermodynamic quantities of the system— in particular the entropy— are 
well defined at these instants. During the processes 2 — 3 and 4 — )■ 1, the rapid change 
(compared with time scale, ~ E/h) in the energy level structure of the system results in 
quantum nonadiabatic phenomenon. We emphasize that in any quantum thermodynamic 
process in which the low dissipation exists the system relaxation is assumed to be fast com- 
pared to the time of the process Sj. During a quantum isothermal (adiabatic) process, the 
entropy production caused by weak dissipation can be written as S/i/ti2 or Sc/t34 (Sa/t23 or 
Sa/i4i), since the reversible regime is approached in the infinite time durations. Thus, the 
real heat exchanges Qh and Qc are = T^AS — Th^ and Qc = T^AS + Tc^ + Tc{^ + f^), 
with AS = S{X2) - S{Xi) = \ [S{X4) - 5(X3)]| - (|f + gf). According to the ffist law of 
the thermodynamics, the work W per cycle is determined by 

W = Q,-Q, = {n- T,)AS - - - T,(^ + ^). (2) 

Let v(t) and r be the speed of the change of X and the cycle period, respectively. Then we 
have 

Xo = (X2-Xi) + (X3-X2) + (X3-X4) + (X4-Xi) 

= 2(X3 - Xi) = / v{t)dt = VT, (3) 
Jo 



where v is the average speed of the change of X. The displacement of X after a single cycle 
is zero and thus X is a state variable, though the total change per cycle Xq is not equal to 
zero. The times spent on two isothermal processes can be expressed as = {X2 — Xi)/v 
and = (X3 — X4)/v, respectively, while the times of two adiabats arc = (^3 — X2)/v 
and = [X^^ — Xi)/v. Therefore, the power output P = W/t and the efficiency rj = W/Qh 
are 

^' r/^ A ^ vThEh vTcEc vTcTja vTcEa n 



2{Xs-X,)'' (X3-X4) (X3-X2) (X4-X1) 

and 



_ l-^ft J-c)^'^ (Xa^Xi) (X3~X4) (X3-X2) (X4-X1) . . 

-^h^^ (X2-X1) 

respectively. Here A^" is a state variable determined only by the initial and final states of 
the isothermal process. The generahzed coordinates Xi with i — 1,2,3,4: , corresponding to 
the system volume V in the classical thermodynamics, are state variables and independent 
of the detailed protocols. To specify the time allocation at maximum power output, the 
values of Xi as well as the average speed v should be optimized. We will the optimize power 
output P over the average speed v and the variables Xi to obtain the time allocation during 
a cycle and thus to determine the corresponding efficiency. We will assume, for simplicity, 
that the initial value of the external parameter is a constant, i.e, Xi — X^. The maximum 
power is therefore found by setting the derivatives of P with respect to the average speed v 
and Xi with i = 2, 3, 4 equal to zero. 

The maximization conditions -^\xi=x^ = and ^\v=vm — gi'^^ the physical solution. 
The value of X^ is determined by the following equation 

(Th - T,)AS = VmTc^ai + ) + + v^TAXl^ - X?) 

\ n cj ^ '^m c a\ vm vm vm' vm vm m c\ ^ 1/ 



^ ^(X^ - X^f {X^ - X^f^ X^ - ^0' n 



where 



(7) 



\-l-h^h^3 — J-c^a^l\ 


+ (^3"^ - 


■ XDy/TcThTio^h 




\Th^h — Tc^a\ 





and 



^ l^a-^3 ~^c^ll + l-^3 --^ijv^a^c ^g-j 



From Eqs. (|6]) and ([7]), we find that the optimal value of with fixed value of is 
independent of the value of the state variable AS*, as expected. Substitution of Eqs. ([7]), 
dH]) and into Eq. leads to the fundamental optimal relationship between X^ and Xf 
at maximum power output. Under the assumption that the value of Xi is fixed at the start 
of the engine cycle, Eq. (jSj) can be done numerically for given values of entropy production 
Sa, Sc and S/j along the specific processes and of temperatures Th and Tc. Once we have 
obtained the optimal relationship between X^ and X™ at maximum power output P, we 
can then determine the optimal values of X^ with i = 2,4, and the average speed Vm by 
Eqs. dZD, dH]), and ©. 

Substituting Eq. ([7]) into Eq. (E]), we find the expression for the efficiency at maximum 
power as follows: 

1 - ^ 

Eq. ( ITOl) together with Eqs. (E]), ([7]), (E]), and ([9]), as one of our main results, conveys the 
following physical features: 

(i) The nonadiabatic dissipation is neglected, i.e., — )• 0. In such a case, the limits 
^ — 7- and ^ — )■ oo, lead to the result that the efficiency r/m at the maximum power 
approaches to the upper bound ?7+ = and to the lower bound r]_ = respectively. 
That is, the efficiency 1]^ at the maximum power satisfies the following condition: 

^^V-<Vm<V^^ TT^- (11) 

In Fig. |2] we plot the efficiency ( ITOl) as a function of r]^ comparing r]^^ with the upper 
and lower bounds (ITT]) . The lower and upper bounds, which are reached in the completely 
asymmetric limits ^ — and ^ — )■ cxd, are identical to the corresponding those derived in 
different approaches However, unlike the previous literature in which the times of two 
adiabats are ignored, the times spent on the two adiabats in the quantum Carnot cycle are 
taken into account. If the symmetric dissipation Sc = S/^ and = \hr', i-e., the time 

allocation to the hot and cold processes at maximum power: 

^12 Th 




tsi y Tc' ^ ^ 

we can recover the CA efficiency rjm = rj^^ = 1 — by using Eq. (ITU]) . Result of 
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consideration of symmetric dissipation agrees with that obtained by optimizing power output 
with respect to the times of the two isothermal processes [s], Q . 

1.0 




FIG. 2: (Color online). Efficiency rjm at maximum power as a function of the Carnot value r]^. 
The upper and lower bounds of the efficiency, ?7+ and ry_ given in Eq. (jlip . are represented by a 
red and a blue solid line with starts and squares, respectively. The CA efficiency ry^^ is denoted 
by a green solid line with 5-pointed stars, while the Carnot efficiency r/^ by a black solid line with 
circles. 



(ii) There is nonadiabatic dissipation, while the dissipation in at least one quantum 
isothermal process is not considered, namely, ^ 0, and S/j — )■ (or Sc — t- 0). From Eq. 
( ITOj) we find in this case 

Vm = V- = \. (13) 

which is independent of the values of both and (or S/i). The lower bound, which is 
found in the case when nonadiabatic dissipation exists but dissipation vanishes in at least 
one isotherm, is particularly interesting. Although it coincides with a reported universal 
lower bound in Refs. js], it is derived in the generalized engine model with nonadiabatic 
phenomenon. Physically, even in the case when the two isothermal processes are reversible, 
inclusion of an arbitrary low dissipation in the adiabats reduces the efficiency at maximum 
power output to half the Carnot value 

(iii) Dissipations in four quantum thermodynamic processes are equal, i.e., Sc/S/i = 1 



and = 1. hetR^ = {X^ - X^){- 



+ 



), in the limits Rx ^ 



and Rx ^ oo, the efficiency rjm corresponding to maximum power output converges to the 
upper bound 77+ = ^^'^ and to the lower bound r]_ = respectively. Here the lower and 
upper bounds are equal to the corresponding those in previous studies, but extended to the 
irreversible QCEs in which the time spent on two adiabats and nonadiabatic dissipation are 
considered. According to Eq. ( ITOl) . the CA efficiency r]^^ = 1 — is achieved when the 
times spent on the four quantum thermodynamic processes are distributed in such a way 
that 

ti2{r - tu) ^ 1% ^^^^ 

where ti2 = {X^ - XD/v^Mz = {X^ - X^)/v^,tu = {X^ - X^)/vm, and Ui = (Xf - 
X^)/u^. 



III. CONCLUSIONS 



In conclusion, we have determined efficiency at maximum power for a QCE engine per- 
forming finite time cycles. To correctly describe the irreversible QCEs the times spent on 
two adiabats and nonadiabatic phenomenon have been taken into account. In the limits of 
extremely asymmetric dissipation — ^ and ^ — )■ 00, with ^ — )■ 0), the efficiency at 
maximum power output converges to an upper and a lower bound, coinciding with the re- 
sult obtained previously in different approaches. When dissipation in any isothermal process 
vanishes, the efficiency at maximum power output is equal to the lower bound For the 
QCE with the symmetric dissipation but without nonadiabatic dissipation — )■ 1, while 
^ — )■ 0), we have derived CA efficiency at maximum power output, only provided that the 
ratio of the times of contact with two heat reservoirs satisfies the relation given as in Eq. 
( !T2l) . In the case of ^ — ?■ 1 and ^ — )■ 1, we have also recovered CA efficiency at maximum 
power output, if the time allocations of four processes fulffil the condition in Eq. (fT^ . 
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